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Self-associates of nucleic acid components (stacking trimers and 
tetramers of the base pairs of nucleic acids) and short fragments 
of nucleic acids are nanoparticles (linear sizes of these particles 
are more than 10 A). Modern quantum-mechanical methods and 
softwares allow one to perform ab initio calculations of the sys- 
tems consisting of 150-200 atoms with enough large basis sets (for 
example, 6-31G*). The aim of this work is to reveal the peculiar- 
ities of molecular and electronic structures, as well as the energy 
features of nanoparticles of nucleic acid components. 

We had carried out ab initio calculations of the molecular struc- 
ture and interactions in the stacking dimer, trimer, and tetramer 
of nucleic base pairs and in the stacking (TpG)(ApC) dimer and 
(TpGpC) (ApCpG) trimer of nucleotides, which are small DNA 
fragments. 

The performed calculations of molecular structures of dimers 
and trimers of nucleotide pairs showed that the interplanar dis- 
tance in the structures studied is equal to 3.2 A on average, and 
the helical angle in a trimer is approximately equal to 30°. The 
distance between phosphor atoms in neighboring chains is 13.1 A. 
For dimers and trimers under study, we calculated the horizontal 
interaction energies. 

The analysis of interplanar distances and angles between nu- 
cleic bases and their pairs in the calculated short oligomers of 
nucleic acid base pairs (stacking dimer, trimer, and tetramer) 
has been carried out. Studies of interactions in the calculated 
short oligomers showed a considerable role of the cross interaction 
in the stabilization of the structures. The contribution of cross 
interactions to the horizontal interactions grows with the length 
of an oligomer. Nanoparticle components get electric charges in 
nanoparticles. Longwave low-intensity bands can appear in the 
electron spectra of nanoparticles. 



1. Introduction 

Self-associates of nucleic acids (stacking trimers and 
tetramers of nucleic acid base pairs) and short NA frag- 
ments are nanoparticles ((linear sizes of these parti- 
cles are more than 10 A). Modern quantum-mechanical 
methods and programs allow one to perform the ab ini- 
tio calculations of systems consisting of 150-200 atoms 
with a rather large basis (6-31G*, for example) [1-3]. 



Earlier, the ab initio calculations were carried out for 
nucleic acid base pairs and for dimers of nucleic acid base 
pairs [4-10]. But the sizes of these structures are less 
than 10 ran. and these structures cannot be considered 
as nanoparticles. Work [11] presents the results of calcu- 
lations for interaction energies of an oligomer including 
7 pairs of nucleic acid bases. But, in those calculations, 
the oligomer geometry was not optimized, being taken 
from NMR data. 

Some time ago, we published the results of calculations 
for the molecular structure of a stacking trimer of nucleic 
acid base pairs [12]. 

The aim of the present work is to reveal the peculiar- 
ities of molecular and electronic structures and energy 
characteristics of nanoparcticles of nucleic acid compo- 
nents. 

2. Subjects and Methods of Calculations 

We performed ab initio calculations of the molecular and 
electronic structures of short DNA fragments - stacking 
(TpG)(ApC) dimer and (TpGpC) (ApCpG) trimer of 
nucleotide pairs with four Na + ions and without Na + , 
as well as stacking dimer, trimer, and tetramer of nucleic 
acid base pairs. 

Geometries of the stacking (AC)(TG) dimer of base 
pairs, (CAC)(GTG) trimer, (ACAC)(TGTG) tetramer, 
and DNA short fragments were optimized by M06-2x 
functional of the DFT method [13]. As was shown in 
[14], this functional describes well the stacking interac- 
tion of nucleic acid components. Then, the geometries of 
dimers and trimers of the base pairs were re-optimized at 
MP2 level of theory. 6-31/G* basis set was used for the 
oligomer geometry optimization. The analysis of bond 
lengths in a dimer calculated by the above methods re- 
vealed a difference of 0.1% in the values of bond lengths. 

The analysis of the vibrational spectra of short 
oligomers showed the absence of imaginary frequencies in 
the calculated spectra of dimers, trimers, and tetramers. 
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TpA...GpC 



Tp...Ap 
Gp...Cp 
C...G + 4Na+ 

Fig. 1. Optimized structures of stacking dimer and trimer of 
nucleotides 

The interaction energies in these structures were deter- 
mined with the typical scheme described in [4]. Accord- 
ing to this scheme, the interaction energy in a complex 
is equal to the difference between the total energy of 
the complex and the sum of total energies of monomers. 
These calculations were performed at MP2 level of the- 
ory, BSSE being taken into account and 6-31+G* basis 
set being applied. The charges on atoms were calculated 
by the Merz-Kolman method [15]. The calculation of 
excited state energies had been carried out by the TD 
DFT method. Here, we used NWChem [1], Gamess [2], 
and Gaussian 03 [3] programs. The visualization of the 
results of calculations is performed by the Chem Craft 
Program [16]. 

3. Results and Discussion 

3. 1 . Molecular structure of DNA fragments 

When considering the calculation results for molecu- 
lar structures of a dimer and a trimer of nucleotide 
pairs (Fig. 1), it should be noted that the calcu- 



lated (ApCpG) trimer is a codon and is able to code 
tryptophan aminoacid [17]. Calculations of the molec- 
ular structure revealed (Table 1) that, in the struc- 
tures studied, the average value of interplanar spac- 
ing is 3.2 A, the helicity angle in a trimer is 30°, 
and distances between phosphorus atoms in neighbor 
chains are 13.1 A, respectively. Intervals between ex- 
treme hydrogen atoms in nucleotide pairs and lengths 
of hydrogen bonds are 15.8 and 1.9 A on the aver- 
age, respectively. Comparison of the calculated cova- 
lent bond lengths in a dimer and a trimer with the X- 
ray structure analysis data revealed their difference of 
about 2%. We calculated the horizontal interactions 
in the dimer and the trimer under study. These ener- 
gies are 40 and 69 kcal/mol, respectively. These ener- 
gies are larger than the total energy of hydrogen bonds, 
and this is related to the contribution of the cross in- 
teractions to the energy of horizontal interactions in 
oligomers of nucleic acid components (as will be shown 
below). 



3.2. Molecular structures of nucleic base pair 
oligomers 

Figure 2 presents the optimized geometries calculated 
for the stacking (AC)(TG) dimer, (CAC)(GTG) trimer, 
and (ACAC)(TGTG) tetramer of nucleic acid base pairs. 

Calculations of the stacking dimer, trimer, and 
tetramer revealed the parameters of structures presented 
in Table 2. 

The interplanar distance in the structures is 3.2 A on 
the average. 

The value of helical twist in the structures lies in the 
— 10° interval (at first, 40° angle was specified). 

The value of propeller twist between bases of the stud- 
ied structures is in the 3-20° interval. The value of pro- 
peller twist in base pairs of the dimer is larger than that 
in the trimer. 

Table 1. Geometry Parameters of Nucleic Acid 



Parameters 


Dimer 


Trimer 


Interplanar distance, A 


3.2 


3.2 


Helical angle, degrees 


10 


30 


Distance between P atoms in 






adjacent chains, A 


13.1 


13.2 


Distance between Na+ ions in 






adjacent chains, A 




15.8 


Horizontal interaction energies, 






kcal / mol 


-40.2 


-69.1* 


* This calculation is performed for a trimer without Na+ ions 



724 



ISSN 2071-0186. Ukr. J. Phys. 2012. Vol. 57, No. 7 



MOLECULAR STRUCTURE AND INTERACTIONS 



Base pairs in the structures are nonplanar. The value 
of dihedral angle between base planes changes within 
160-180°. 

Middle pairs in the trimer and in the tetramer are 
more planar (the dihedral angle and the propeller one 
are equal to 0) than outside ones. On the whole, pairs in 
the trimer are more planar than those in the dimer. This 
is an important peculiarity of the molecular structure 
calculated for the trimer and the tetramer. 

3.3. Interactions in stacking oligomers of base 
pairs 

We performed calculations of the vertical and horizontal 
interactions in the oligomers under study. The vertical 
interaction energies between AT and CG pairs are equal 
to —12.02 kcal/mol in a dimer and —10.98 kcal/mol (be- 
tween high and middle pairs) or —11.0 kcal/mol (be- 
tween middle and down pairs) in a trimer, respectively. 
These data correlate with results of work [2], in which 
the parameters of the vertical interactions in stacking 
dimers of base pairs were calculated by the complete ba- 
sis set extrapolation. 

It is of interest to evaluate the energies of horizontal 
interactions in the studied structures. Table 3 presents 
the calculation results for the energies of horizontal inter- 
actions in the oligomers and the sums of hydrogen bond 
energies in AT and GC pairs forming these oligomers. 
Our calculations of the energies of hydrogen bonds in 
isolated pairs of AT and GC gave —13.83 kcal/mol for 
AT and —26.98 kcal/mol for GC pairs. As is seen from 
Table 3, the energies of horizontal interactions exceed 
the arithmetic sum of the energies of hydrogen bonds in 
isolated AT and GC pairs. These overridings are equal 



Table 2. Structural parameters of calculated stacking 
dimer, trimer, and tetramer 



Structural 
parameter 


(AC)(TG) 


(CAC)(GTG) 


(ACAC)(TGTG) 


Helical twist, 


0° 




0° 


CG1-CG3 


5° 


AT1-CG2 


degree 






10° 


CG1-AT2 


2° 


CG2-AT3 








10° 


AT2-CG3 


6° 


AT3-CG4 


Propeller twist, 


15° 


AT 


5° 


CGI 


2° 


ATI 


degree 


20° 


GC 


0° 


AT 


11° 


CG2 








3° 


CG3 


0° 


AT3 












8° 


CG4 


Dihedral angle 


175° 


AT 


178° 


CGI 


161° 


ATI 


between bases, 


176° 


GC 


180° 


AT 


176° 


CG2 


degree 






177° 


CG3 


179° 


AT3 



-170° CG4 




(CAC)(GTG) 




(CACA(GTGT)) 

Fig. 2. Optimized structures of the stacking (AC) (TG) dimer, 
(CAC) (GTG) trimer, and (ACAC) (TGTG) tetramer of nucleic 
base pairs 

to —2.7 kcal/mol, —9.39 kcal/mol, and —13.1 kcal/mol 
for the dimer, trimer, and tetramer, respectively. 

To explain the contradiction revealed between the 
energies of horizontal interactions and the sum of the 
energies of hydrogen bonds in base pairs forming the 
oligomers, we used the idea of cross interactions [18] and 
calculated the energies of cross interactions in the sys- 
tems under study. Cross interactions are induced by 

Table 3. Energies of horizontal interactions (kcal/mol) 
in the stacking (AC)(TG) dimer, (CAC) (GTG) trimer, 
and (ACAC) (TGTG) tetramer calculated by various 
methods (6-31+G* basis set) and the sums of the energies 
of H-bonds (kcal/mol) in isolated base pairs 



Oligomer/ 
Method 


MP2//MP2 


MP2//M06-2X 


Sums of H-bond 
energies 


Dimer 


-42.62 


-43.51 


-40.81 


Trimer 




-77.18 


-67.79 


Tetramer 




-94.72 


-81.62 
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Fig. 3. Calculated scheme of interactions in the trimer. Assig- 
nations: H-b - hydrogen bonds, St -stacking interactions, Cr - 
cross-interactions 

Coulomb, dispersion, and inductive interactions and in- 
clude interactions of every present base with not only 
the nearest base, but with other ones in a stack or in a 
helix (Figure 3). 

So, the horizontal interaction energies must include H- 
bond energies and cross interaction energies. The verti- 
cal interaction energies must include stacking interaction 
and cross interaction energies. 

3-4- Dimer 

Table 4 presents the calculation data on the energy of 
horizontal interaction in the dimer of base pairs, energies 
of hydrogen bonds in base pairs in the dimer configura- 
tion, energies of the vertical interaction between AT and 
CG base pairs, energies of stacking interactions in the 
dimers of bases, being a part of the base pair dimer, as 
well as the energies of interactions between bases, being 
a part of different pairs (cross energies). In all these cal- 
culations, the geometry of bases was taken into account, 
which they possessed in the base pair dimer composition. 
As stated above, the dimer geometry was optimized at 
MP2 level of theory. 

The following designations are accepted in Table 4. 
The second column (lines 1 and 4) presents the schemes 
of relative positions of nucleic acid bases - adenine (A) , 
thymine (T), guanine (G), and cytosine (C) in the cal- 
culated stacking dimer of nucleic acid base pairs. This 
column also shows the schemes of relative positions of 
nucleic bases in (A)(T), (G)(C) base pairs, in stacking 
dimers of (A)(C) and (T)(G) pairs, and in cross pairs 
(A)(G) and (T)(C). Horizontal arrows (line 1) mark the 
horizontal interaction breakdown in the stacking dimer 
of base pairs. Horizontal arrows in lines 2 and 3 point 
to the opening of hydrogen bonds in (A)(T) and (G)(C) 



base pairs. Underlining in schemes (lines 4, 5, 6) shows 
the stacking interaction failure in the dimers. The third 
column presents the interaction energies calculated for 
the structures shown in column 2. Column 4 gives the 
summation results for hydrogen bond energies and cross 

Table 4. Energies of interactions between compo- 
nents inside the stacking (AC)(TG) dimer (MP2//MP2 
method, 6-31+G* basis set) 



No. 


Schemes" of structures 
calculated 


Interaction energies, 
kcal / mol 


Sums 9 , 
kcal / mol 


1 


(A)^(T) 6 


-42.46* 


-42.86 




(C)«(G) 




2 


(A)(T) C 


-12.80 




(...)(...) 




3 


(-)(-) c 


-26.31 




(C)«->(G) 




1 


(A)(T)<* 


-12.02" 


-12.39 




(C)(G) 




5 


(A)(.„r 


-3.88. 




(C)(...) 




6 


(...)(T)« 


-4.77 




(••0(G) 




7 


(A)(...)/ 


-3.73. 




(-O(G) 




8 


(-•OCT)' 


-0.01 




(C)(...) 





a Scheme of the (AC)(TG) stacking dimer in this Table is ( A ^ T ) 

(C)(G) 

(A) -o- (T) 

b Scheme means the breaking of the horizontal inte- 

(C) «• (G) 

raction in the (AC)(TG) stacking dimer. 

(A) (T) ( )( ) 

c Schemes and v""/v-v show the breaking of the 

(...)(...) (C)^(G) 
H-bond interaction in AT and CG pairs in the (AC)(TG) 
stacking dimer. Geometries of these AT 

and GC pairs are identical to ones in the studied dimers. 
(A)(T) 

d Scheme shows the breaking of the stacking interaction 

(C)(G) 

between AT and GC pairs in (AC)(TG) stacking dimer. 

(A)(...) (...)(T) 
e Schemes and show the breaking of the sta- 

(C)(...) (...)(G) 
eking interactions in the stacking AC and TG dimers as a 
compo- nent of the (AC)(TG) stacking dimer of base 
pairs. The geomet- 
ries of these AC and TG dimers are identical to ones in the 
(AC)(TG) dimer. 

(A)(...) (...)(T) 
f Schemes and show the breaking of cross in- 

(-O(G) (C)(...) 
teractions in the (AC)(TG) stacking dimer between A and G and 
between T and C, respectively. 

9 Sums of all possible interactions (H-bonds, stacking and cross 
interactions) of components are shown. 
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interaction energies (line 1) and for stacking and cross 
interaction energies (line 4) in the calculated base pair 
dimer. It should be noted that the energy of hydrogen 
bonds in AT and GC pairs being a part of the dimer is 
lower than that of hydrogen bonds in isolated AT and 
GC pairs, which is related to the different geometries 
of these pairs in the dimer composition and in the free 
state. 

As Table 4 shows (column 3), the energy of horizontal 
interaction in the dimer is by 3.35 kcal/mol higher than 
the arithmetic sum of the energies of hydrogen bonds 
in AT and CG base pairs in the dimer. The energy of 
vertical interaction in the dimer between AT and CG 
pairs is more than the arithmetic sum of the interaction 
energies in the dimers of AC and TG bases, being a part 
of the stacking base pair dimer. 

As is noted above, this difference is probably caused 
by the cross interactions. Really, the calculations of 
the interaction energies between bases lying in differ- 
ent planes (cross energies) result in the values of —3.73 
kcal/mol and —0.01 kcal/mol for AG and TC pairs, re- 
spectively. Taking these interactions into account allows 
one to reach a rather good agreement between the hori- 
zontal and vertical interaction energies and the sums of 
energies between monomers in the stacking dimer (most 
right column in Table 4). Differences between the total 
energy calculated for horizontal interactions in the dimer 
and the arithmetic sum of hydrogen bonds in base pairs 
and of cross interactions are within 1%, similar to the 
difference being about 3% for the vertical interaction. 

So, the presence of the cross interaction explains a 
rise of the horizontal and vertical interaction energies in 
the stacking dimer, in comparison with the sum of the 
energies of hydrogen bonds in base pairs and the sum 
of the energies of stacking interactions in the dimers of 
bases. 

3.5. Interactions in the trimer and the tetramer 

The same situation is observed both for the trimer and 
the tetramer. Table 5 presents the energies of horizontal 
interaction in the trimer of base pairs, energies of hydro- 
gen bonds in base pairs in the trimer, energies of vertical 
interactions between AT and CG base pairs, and ener- 
gies of stacking interaction in the dimers of bases, being 
a part of the trimer, as well as the energies of interactions 
between bases forming different pairs (cross interaction 
energies). In all these calculations, the geometry of bases 
was used, which they have in the trimer structure. In 
addition, Table 5 presents the results of summation of 
hydrogen bond energies (rightmost column) and of the 



energies of AG and TC cross interactions in the trimer 
(above) and the results of summation of stacking interac- 
tion energies in the dimers of AC and TG bases forming 
the studied stacking trimer of (CAC)(GTG) base pairs 
and the energies of AG and TC cross interactions. 

It should be noted that the energies of hydrogen bonds 
in AT and GC pairs in the trimer differs both from those 
of hydrogen bonds in the dimer and in isolated AT and 
GC pairs, which is determined by a difference in the ge- 
ometries of these pairs in the trimer, dimer, and in the 
free state (Table 2). As is seen from Table 5, the calcu- 
lations revealed that the energy of horizontal interaction 
in the stacking trimer exceeds the arithmetic sum of the 
hydrogen bond energies in isolated AT and GC pairs. 
The excess makes up —10.84 kcal/mol. Energies of the 
vertical interactions in the trimer between AT and CG 
pairs are larger than the arithmetic sums of interaction 
energies in the dimers of bases - AC and TG being a 
part of the stacking trimer composition. 

The calculations of interaction energies of the bases 
placed in various planes (cross interaction energies) re- 
sulted in —4.39 kcal/mol and —4.47 kcal/mol for GA 
and AG pairs and —0.63 kcal/mol for CT and TC pairs, 
respectively. There is some difference of interaction en- 
ergies for reciprocal GA and AG pairs. Table 5 also 
presents the energies of cross interactions between bases 
(adenine and guanine) in the GAG triad and cytosine 
and thymine in the CTC triad of the stacking trimer 
under study. Just taking these cross interactions into 
account allows us to reach a very good agreement be- 
tween the total energy of horizontal interactions and the 
sum of the energies of interactions between monomers 
(nucleic acid bases) in the trimer. 

Table 5 shows the energies of vertical interactions be- 
tween the upper GC and middle AT pairs and between 
middle AT and down GC pairs, which are only by about 
0.7 kcal/mol lower than the summation results for stack- 
ing and cross interactions of components in these struc- 
tures. Thus, as in the dimer case, the consideration of 
the cross interaction explains a rise of the energy of hor- 
izontal and vertical interactions in the base pair trimer 
in comparison with the sum of the energies of hydro- 
gen bonds in base pairs and with the sum of stacking 
interactions in dimers of bases, being the part of the 
trimer. 

Let us consider the interaction energies for the 
tetramer. Table 1 from Supplement presents the cal- 
culation results (column 3) for the total energy of hori- 
zontal interaction in the tetramer, energies of hydrogen 
bonds in base pairs being a part of the tetramer, and 
energies of paired interactions (stacking and cross in- 
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Table 5. Energies of interactions between components inside the (CAC)(GTG) trimer (MP2/ /M06-2X method, 
6-31+G* basis set) 



No. 


Schemes of structures 


Interaction energies, 


Sums b , 


No. 


Schemes of structures 


Interaction energies, 


Sums b , 




calculated a 


kcal / mol 


kcal / mol 




calculated a 


kcal / mol 


kcal / mol 



(C)4 

(A)< 
(C)+ 
(C)+ 



(A). 



KG) 
>(T) 
KG) 

KG) 
(...)(...) 
(...)(...) 
(...)(...) 
(...)(...) 
(C)<-f(G) 
(...)(...) 
KT) 

(...)(...) 

(C)(G) 
(A)(T) 
(C)(G) 
(C)(G) 
(A)(T) 
(C)(G) 
(C)(...) 
(A)(...) 
(...)(...) 
(...)(...) 
(A)(...) 
(C)(...) 



-77.18 



10 d 



-26.61 

10.0 
-13.11 



-11.00 



-10.98 



-2.51 



-2.43 



-76.46 c 
-77.16 d 

0.8170 



0.87273 



-11.70 



-11.72 



10 



11 



12 



13 



14 



15 



16 



(...)(G) 
(■■■)(T) 
(•••)(•••) 
(•••)(•••) 
(■■■)(T) 
(...)(G) 

(•••)(•••) 
(A)(...) 
(...)(G) 
(...)(G) 
(A)(...) 
(•••)(•••) 
(■■■)(G) e 
(A)(...) 
(...)(G) 
(C)(...)/ 
(■■■)(T) 
(C)(...) 
(•••)(•••) 
(■■■)(T) 
(C)(...) 
(C)(„.) 
(■■■)(T) 
(•••)(•••) 



-4.17 



-4.19 



-4.47 



-4.39 



-9.38 



-1.44 



-0.63 



-0.63 



(C) o (G) 

"Scheme of the (CAC)(GTG) stacking trimer in this table is (A) (T). Other designations in this Table are similar to those in 

(C)++(G) 

Table 4. Horizontal interaction in this Table is between (CAC) and (GTG) stacking trimers of nucleic bases. Stacking interactions are 

between the upper CG pair and (AC)(GT) components and between (CA)(GT) and the down CG pair. 

'Sums of all possible interactions (H-bonds, stacking and cross interactions) of the components. 

c Sum of H-bond interactions with using separate GA, AG, CT and TC (lines 11, 12, 14, 15) cross interactions. 

d Sum of H-bond interactions together with GAG (linel3) and CTC (line 16) cross interactions. 



teractions) in the tetramer. As in the cases of dimer 
and trimer, the energies of horizontal and vertical in- 
teractions exceed the arithmetic sum of the energies of 
hydrogen bonds and stacking interactions in pairs and 
the dimers of bases being in the tetramer. Nevertheless, 
the consideration of the cross interactions in cross pairs 
and triads of bases present in the tetramer allows us to 
"compensate" this difference. 

It is worth noting the variability of the energies of 
cross interactions in AG and TC pairs in the tetramer. 
The energies of cross interactions change from —4.12 
kcal/mol to —4.70 kcal/mol and from —0.42 kcal/mol 
to —0.81 kcal/mol for AG and CT pairs, respectively. 
This variability is related to the difference in the spatial 



structures of complementary AT and GC base pairs (in 
particular, with the difference in dihedral and propeller 
angles, which include AG and CT "cross pairs" (see Ta- 
ble 2)). 

In addition, we calculated the energies necessary for 
the "opening" of end AT and GC base pairs, being a part 
of the tetramer. The values of energies were —23.36 
kcal/mol and —37.49 kcal/mol for AT and GC pairs, 
respectively. The values of energies calculated directly 
for end pair openings are close to those of the sum of 
the energies of hydrogen bonds in base pairs, as well as 
the energies of cross and stacking interactions between 
the nearest bases of nucleic pairs in the tetramer. The 
calculations of the opening energy of the outer GC pair in 
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Fig. 4. Molecular orbitals (HOMO and LUMO+6), between which 



Table 6. Total charges of base stacks in the 
(AC)(TG) dimer, (CAC)(GTG) trimer, and (AC AC) 
(TGTG) tetramer 



Stacks of nucleobases 
in oligomers 


Total charges in stacks of 
nucleobases 


(AC) 


0.115 


(TG) 


-0.115 


(CAC) 


0.229 


(GTG) 


-0.229 


(ACAC) 


0.243 


(TGTG) 


-0.243 



the tetramer and the trimer gave similar results: —37.49 
kcal/mol and —37.98 kcal/mol for the tetramer and the 
trimer, respectively. 

The analysis of Tables 3-5 and Table 1 in Supple- 
ment allowed us to reveal increasing the relative contri- 
bution of cross interactions to the total horizontal in- 
teractions in the structures under study. This contribu- 
tion increases from 7.9% for the dimer to 18.4% for the 
tetramer. 

3. 6. Charges 

The analysis of the electronic structure of oligomers (in 
particular, the analysis of the electronic density distri- 
bution) is of special interest. The calculations of the 
charges on atoms (Table 6) showed that the sum of 
charges in stacks of (AC)and(TG), (CAC) and (GTG), 
(ACAC) and (TGTG) bases forming the oligomers un- 
der study is not equal to and comes to 0.115, 0.229, 
and 0.243 for the dimer, trimer, and tetramer, respec- 
tively. It is a manifestation of the charge transfer in the 
horizontal direction in the short stacking oligomers. 

3.7. Excited states 

The analysis of excited states in the structures is of sig- 
nificant interest as well. The calculations of the excited 




LUMO+6 

the first electronic transition occurs in the dimer 



Table 7. Parameters of dimer (left) and trimer (right) 
excited states calculated by the TD DFT method on MP2 
level of theory calculated geometries 



Excited 


Excited 


Wave 


Oscil. 


Excited 


Wave 


Oscil. 


state 


state 


length, 


strength 


state 


length, 


strength 




energies, 


nm 




energies, 


nm 






eV, dimer 






eV, trimer 






1 


4.70 


263.4 


0.0040 


4.60 


269.1 


0.0030 


2 


4.86 


254.9 


0.012 


4.64 


267.1 


0.0034 


3 


4.91 


252.6 


0.038 


4.69 


264.2 


0.0033 


4 


4.94 


250.1 


0.016 


4.75 


260.7. 


0.0479 


5 


5.05 


245.3 


0.001 


4.76 


260.3 


0.0004 


6 


5.12 


242.3 


0.008 


4.81 


258. 


0.0401 


7 


5.14 


241.1 


0.050 


4.82 


257.3 


0.0051 


8 


5.17 


239.7 


0.282 


4.82 


257.0 


0.0132 


9 


5.24 


236.3 


0.101 


4.95 


250.0 


0.001 


10 


5.36 


231.2 


0.243 


4.99 


248.21 


0.0558 



state energies in the stacking dimer and trimer revealed 
a possibility for a low-intensity band to exist in the long- 
wave region of the spectrum (Table 7). This corresponds 
to our previous suppositions [19] on possible charge- 
resonance interactions in crystals and dimers of nucleic 
acid bases, as well as to the data of work [20]. The 
analysis of molecular orbitals in the dimer showed (Fig. 
4) that this low-frequency transition has to take place 
between the molecular orbital localized on guanine and 
the molecular orbital localized on the neighboring pair 
of bases (thymine and adenine). Thus, we deal with the 
electronic transition between two systems of molecules 
being placed one over another. Moreover, in the trimer 
case, we can observe several low-intensity electronic tran- 
sitions in the low-energy region of the spectrum. 

4. Conclusions 

1. The ab initio calculations of nanoparticles of nucleic 
acids components - a dimer and a trimer of nucleotides 
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Table 8. Energies of interactions between components 
inside a tetramer (MP2//M06-2X computation method) 



No. 




Tn1"pT';}p1"ini'i 


No. 


1 1*3 1 1 1 o 






^trncMiTP^ 


pnprci 




structures 


PnPT'0"lPS 




Cctlculfi-trGCl 


kcal / mol 




Ccilculcttcd 


kcal / mol 


x 


K-rt-)^ * V s -) 


—94.72 


9 




0.42 










(cm i 






(A)< — ¥ (T) 






( )( ) 






(C)< y (G) 






( )( ) 




2 


( )( ) 


-26.79 


10 


f XT) 


-0.81 




(C,}4 y <G) 












V. — / \ — / 






( )(T) 












(...)(...) 

\"'/ \ / 




3 


(...)(...) 


26.44 


11 


(...)(...) 


-0.76 




(■■■)(■■■) 

\"'/\"'/ 






(C)(...l 






(...)(...) 












(C)< — >■ (G) 






(C)(...) 




4 


(A)<— y (T) 


-12.03 


12 


(A)(...) 


-2.59 




(•••)(•••) 






(C)(...) 






(•••)(•••) 






(■■■)(■■■) 






(•••)(•••) 






(...)(...) 




5 


(•••)(•••) 


-12.04 


13 


(■■■)(■■■) 


-2.446 




(•••)(•••) 






(■■■)(■■■) 






(A)-i — > (T) 






(A)(...) 






(•••)(•••) 






(C)(...) 




6 


<A)(...) 


-4.70 


14 


(...)(...) 


-3.64 




(...)(G) 






(...)(...) 






(•••)(•••) 






(-)(T) 






(...)(...) 






(...)(G) 




7 


(•••)(•••) 


-4.35 


15 


(-)(T) 


-2.55 




(•••)(•••) 






(,..)(G) 






(A)(...) 






(•••)(•••) 






(...)(G) 






(...)(...) 




8 


(•••)(•••) 


-4.12 










(...)(G) 












(A)(...) 












(•••)(•••) 











Assignations are the same as in Tables 4 and 5. 



and a stacking trimer and a stacking tetramer of nucleic 
bases have been carried out for the first time. 

2. The calculations performed revealed peculiarities of 
the structure and the interaction energy of nanoparticles: 

- Cross interactions play a considerable role in the sta- 
bilization of the structures. The contribution of cross 
interactions into horizontal interactions grows with the 
lengthening of the oligomer. 

- Internal pairs of bases in a trimer and a tetramer of 
base pairs are more planar in comparison with external 
ones. 

- Nanoparticle components get electric charges in 
nanoparticles. 



- Low-intensity long-wave bands can appear in the elec- 
tronic spectra of nanoparticles. 
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MOJIEKYJlflPHA CTPYKTYPA TA B3A6MOAII 
y HAHOHACTHHKAX KOMIIOHEHTIB 
HYKJIEIHOBHX KHCJIOT: HEEMIIIPHHHI 
P03PAXYHKH 

W.B. Py6in, Jl.r. Benoyc 
P e 3 io m e 

CaMoacouiaTH KOMnoHeHTiB HyKjiemoBiix khcjiot (cTSKHHr Tpmie- 
pn Ta TeTpaMeTpH ochob HyKjiemoBiix khcjiot) Ta kopotki cppa- 
rivieHTH HyKjieiHOBHx khcjiot e HaHonacTHHKaMH (i'xHi jiiHiimi po3- 
Mipn 6ijn>Hie 10 A). Cynacm KBaHTOBO-MexamnHi mcto/th Ta npo- 



rpaMH JT03BOJI5HOTB 3iJ,iHCHIOBaTH p03paxyHKH CIICT6M, IHO M1CT5ITB 

150—200 aTOMiB, i3 jjpcTaTHBO bcjikkhm 6a3HCOM (HanpHKjia^, 6- 

31G*). MeTOK) po6oTH e BHHBJieHHH OCo6jIHBOCTeft MOJieKyjIHpHO'i 

Ta ejieKTpoHHoi CTpyKTyp, eHepreTHHHHx xapaKTepncTHK EaEOia- 

CTHHOK KOMnOHeHTiB HyKJieiHOBHX KHCJIOT. 

HaMH BHKOHaHO HeeMnipH^Hi po3paxyHKH MOjieKyjrapHOi CTpy- 
KTypH Ta B3aeMOjj;iH b CTeKHHr jxHMepi Ta TpHMepi nap HyKjieoTH- 
jj,iB, a TaKO»c b CTeKHHr jTHMepi, TpHMepi Ta TeTpaMeTpi nap ochob 

HyKJieiHOBHX KHCJIOT. 

ItpoBeflem po3paxyHKH MOjieKyjinpHoi CTpyKTypn jj,HMepy Ta 
Tpnuiepy nap HyKjieoTHjxiB noKasajiH, mp MiscnjionniHHa Bijj,CTaHB 
b CTpyKTypax, mp BHBHajiHCH, b cepejxHBOMy jjppiBHiOBajia 3,2 A, 
KyT cnipajiBHOCTi b TpHMepi - 30°, BijrcTaHB mijk aTOMaMH cpo- 
ccpopy b cycijxmx jiaHLuorax 13,1 A. TaKO»c po3paxoBam eHepri'i 
ropH30HTajiBHHx B3aeMOjxift y nirx /niMepax Ta Tpninepax. 

Mh 3Po6hjih fleTajiBHHit aHani3 Mi>KnjioiHHHHHx BijjpTaHeft i Ky- 
tib Misc ocHOBaMH Ta ix napainH y CTeKHHr jxHMepi, TpHMepi Ta 
TeTpaMeTpi nap ochob HyKJieiHOBHX khcjiot. AHajii3 ropn30HTajiB- 
hhx i BepTHKajiBHHx B3aeMOjj,ijS b po3paxoBaHHx kopotkhx ojiiro- 
Mepax nap ochob HyKJieiHOBHX khcjiot noKa3aB 3HanHy pojiB ne- 
pexpecHHx B3aeMOjj;iH y CTa6ijii3aHi'f biibeghhx CTpy KTy p . Bhccok 
nepexpecHHx B3aeMOjj,iit b ropH30HTajiBm B3aeMOjxi'i 36ijiBHiyeTBCH 
3 noflOBjKGHHHM ojiiroMepy. IIoKa3aHO, mp y CKjiajri HaHonacTHHOK 

lXHi KOMnOHeHTH OTpHMyiOTB ejieKTpHHHHH 3apHJ2. B ejieKTpOHHHX 
CneKTpax HaHOHaCTHHOK MO)KyTB 3'hBHTHCB flOBTOXBHJIBOBi Majio- 
iHTeHCHBHl CMyrn. 
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